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Abstract 

I We investigate the hyperfine sphtting of the heavy baryons in the bound-state 

^ \ approach. We start with an ordinary relativistic Lagrangian which has been extensively 

used to discuss finite mass corrections to the heavy limit predictions. It turns out that 
the dominant contribution arises from terms which do not manifestly break the heavy 
spin symmetry. The actual heavy spin violating terms are uncovered by carefully 
performing a 1/M expansion of this Lagrangian. 
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1 Introduction 



There has been a great deal of recent interest in studying heavy baryons in the bound-state 
picture [0, together with heavy-quark spin symmetry 0. This approach raises many 
fascinating questions which have been explored by several groups P-[TO|. 

These models consist of a chiral Lagrangian for the light flavors and a Lagrangian, £heavy, 



which contains the heavy-meson multiplet H. The simplest choice for the latter is [|TT 

>Chcavy/M = iV^ Tr [hD^h] + id Tr [h^^^^p^H 

where is the four-velocity of the heavy particle and = — iv^ is the covariant chiral 
derivative. Furthermore = (V2) i^d^^"^ ± ^d^^, wherein ^ = exp {icjy/FT,) is the non- 

linear representation of the light pseudoscalar mesons 0. Finally M is the heavy meson mass 
while d is a heavy meson-light meson coupling constant. The light part of the Lagrangian 
allows for a soliton configuration ^c- In the bound state approach the heavy baryon then 
emerges as a heavy meson bound state in the background of ^c- The predictions are very 
simple in the limit where both Nq and M go to infinity. For example, the binding energy 
of the heavy baryon H, |[ is (3/2)(iF'(0) where -F'(O) is the slope of the soliton profile at 
the origin. 

An immediate question is how to estimate what happens when we consider realistic values 
for M. In general this requires the addition of many unknown terms to Eq. (|^). A predictive 
model for finite M corrections may be obtained by constructing a Lagrangian C [0, |1^] of 
a heavy pseudoscalar meson P and a heavy vector meson Q^: 

+ 2iMd (Pp^Q^ - Qf^Pf^P) - id'e^pf,, (D^Qpp^Q, - Qo^ppD^Q,) . (2) 

This reproduces Eq. (|T]) for large M when d' = d and M* = M. We have used D^P = 
{d^ — iv^)P and Q^^ = {jD^Qj^ — D^Q^^. This model in particular allows for manifest 
breaking of the heavy spin symmetry by choosing M* ^ M and/or d! ^ d. The Lagrangian 
(0) represents the starting point for computing physical quantities along the lines of the 
original bound state approach |I|] to strangeness in the Skyrme model p!3| , [1^. This requires 



the solutions to the equations of motion for P and in the soliton background. The 
calculation |10] exhibits sizable corrections for finite M. In addition, recoil effects (finite 
Nq) seem to be very important as well p], When both these effects are taken into 



account it becomes difficult to fit the existing experimental data on the spectrum of the 
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lOH that the inclusion of hght vector mesons 



heavy baryons. It was, however, noticed 
appreciably improves the situation. 

In this note we will resolve an apparent puzzle which arises when calculating the correc- 
tions to the hyperfine splitting using Eq. (□). 



2 An apparent puzzle 

First let us consider the calculation of the hyperfine splitting in the heavy field approach. 
This, of course, arises at first sub- leading order in 1/M and violates the heavy spin symmetry. 
Thus we must add to Eq. (|l]) suitable heavy spin violating terms 0: 

£'h,,,y/M = Tr [Ha^jla,,] - ^{d - d')Ti [Hp;H-i,^,] +■■■ . (3) 

The first term has no derivatives while the second term has one derivative. The hyperfine 
splitting is related to a collective Lagrangian parameter (see section 4 for details) x with a 
proportionality factor of the A-iV mass difference: 

m(E^) - m(SQ) = [m(A) - m{N)] x ■ (4) 

(At present only Sc is well established experimentally.) For Eq. (^) we have 

_M*-M d-d' 
^ ~ 4dF'{0) ^ ~Jd~ ■ ^ ' 

The first term was obtained in Ref. |0] while the second seems to be new. Notice that 
(M* — M) and {d — d') behave as 1/M. These quantities are the same as the ones appearing 
in the ordinary field Lagrangian (|^). It would thus seem that C'^^^^.^^ in Eq. (0) neatly 
summarizes the heavy spin violation in Eq. (Q). 

Now let us consider the calculation of x from Eq. (^ directly based on exact numeri- 
cal solution of the associated coupled differential equations. We content ourselves with the 
graphical presentation of some results^ and relegate the details to a forthcoming publica- 
tion Figure |^ shows x plotted against M for three cases: i) M* = M, d' = d = 0.53, 
ii) M* - M (0.258GeV)VM (a fit to experiment), d' = d = 0.53, iii) M* = M, 
d' — d = (0.0991GeV)/M (an arbitrary choice which sets the coupling constant splitting 
to be 10% at the D meson mass). We immediately notice that x does not vanish when there 

* For the Skyrme model parameters we use the experimental value of i<V and esk = 6.0. This results in 
a profile with F'{0) = f.20GeV. 

t Similar calculations were done in Ref. M but they did not consider the M = A-I*, d ~ d! case. 
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Figure 1: x vs. M computed by numerical integration. Solid line M* = M, 
d! = d; dotted line M* ^ M, d' = d, dashed line M* = M, d' ^ d. 

is no manifest heavy spin violation, i.e., M = M*, d = d' . In fact the dominant part of 
the contribution to x for realistic heavy meson masses is already present in this case. By 
subtracting out this piece we note that the signs of the contributions due to M* ^ M and 
d' ^ d agree with those predicted in Eq. (j^). It is interesting to note that all three curves 
in Fig. |l] fall off as 1/M for M > 10 GeV. But our main task is to understand the source of 
the puzzling non-zero contribution in case i. It is clear that the ordinary field Lagrangian 
(@) must contain heavy spin violating pieces which are not manifest. We will now explore 
this in detail by rewriting Eq. in terms of the "fluctuation field" H and expanding it in 
powers of 1/M. 

3 Expansion of Lagrangian 

Since the effects of M 7^ M* and d ^ d' were taken into account in Eq. (j^) it is sufficient 
to expand Eq. (0) with M* = M and d' = d. To describe the heavy particle moving with 
four-velocity V^, we introduce the factorization 

P = e^^^-^P' , = e^^^-^g^, . (6) 
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P' is the pseudoscalar "fluctuation field" . is not exactly the vector fiuctuation field since 
V ■ Q is not constrained to be zero. We therefore introduce the correct fiuctuation field Q'^ 

by 

Q, = Q',~V,V-Q , (7) 

which shows that V ■ Q' = 0. Substituting Eqs. (|^) and (]^) into the Lagrangian gives, in 
addition to the leading terms of order M, the presently interesting terms of order M": 



+ MV -QV-Q-iM [D^Q'^V -Q-V- QD^Q'^ 
- 2iMd {p'V - pV -Q - V -QV ■ p'P^ + ■ ■ ■ , (8) 

where the three dots stand for terms of order 1/M. In contrast to the massless fields P' and 
Q', V ■ Q is seen to have the large mass M. We thus integrate it out using the equation of 
motion 



(9) 



Substituting Eq. (H) back into Eq. (H) gives 

£(P,g) = {oidei M)+P'D^P^ + Q'^D^Q^-iQ'^F^,{v)q^ 
- 2d (p'V ■ pD^Q^ + D.Q'^V ■ p-p) 

-U^P' [V -pfT^ +■■■ , 



(10) 



where F^y{v) = d^v^, — dyV^ — i[Vfj_,Vy\. In order to extract the heavy spin violating pieces 
it is convenient to rewrite the order Lagrangian in terms of the heavy multiplet field 



H 



'1 — i'^ ■ V) (75P' + 7 ■ Q'). After some algebraic calculation we find 



C{H) 



heavy 

+ d 



-Tr 



HD^H 



+ i-Tr 



2 



D^H^^^,iV-p)H 



- -Tr 
4 



7 ■ DH-i^-i^p^H 



+ d^ 



- -Tr 

4 

It. 

2 



7 ■ DHp^H-i^-i5 + -Tr a^^DaH-faV ■ pi-,(T^yH + h.c. 



H {V ■ pf h] + ^Tr [a^,Ha^, {V ■ pf H 



+ 



fill 



where £heavy is given in Eq. ([1|). At this stage we see that Eq. ( [Tl| ) actually contains pieces 
which are not manifestly invariant under the heavy spin transformations H —>■ SH, H —>■ 
HS^ . These pieces involve two derivatives. 
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4 Hyperfine splitting 



We now sketch the computation of the portion of x in Eq. (^) which results from the "hidden" 
heavy spin violation in Eq. (0) that has been made explicit in Eq. (pJ]). For this purpose 
one needs the collective Lagrangian of the quantum variable A[t) which is obtained after 
substituting 

i{x, t) = A{t)Ux)A\t) , H{x, t) = A{t)H:{x) , (12) 

(where C,c{^) is the classical Skyrme soliton and Hc{x) is the heavy meson bound-state wave 
function) and integrating over d^x. The key dynamical variable is the "angular velocity" fl 
defined by A^A = |t ■ f2. The bound-state wave function may be conveniently presented in 
the rest frame where 

"-(;.:)■ 

with a, I, h representing respectively the iso-spin, light spin and heavy spin bivalent indices. 
We write 

where u{r) is a radial wave function (assumed very sharply peaked near r = for large M) 
and, to leading order in M, the "angular part" of the ground state wave function is 0, |^ 

V'ili = ■^^'iiS2h ■ (15) 

The specific value of the index h results from the choice G3 = G = 1/2 where G is the 
"grand spin" . To next leading order in M the ground state wave function receives a heavy 
spin violating admixture of 

Finally, the hyperfine splitting parameter x is recognized by expanding the collective La- 
grangian [|l|, in powers of O and picking up the linear piece Lcou = (x/2)f23 -|- • ■ ■. Noting 
that the A-nucleon mass difference is given by the moment of inertia, which relates the 



angular velocity to the spin operator |T^, this piece of the Lagrangian yields Eq. (^ after 
canonical quantization of the collective coordinates There are two types of contribution 
to X- The first type, from the heavy spin violating terms proportional to d in Eq. (pTj), 
corresponds to the evaluation of heavy spin violating operators in the ground state (|15]). 
The second type corresponds to the evaluation of heavy spin conserving operators in the 
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ground state which includes an admixture of Eq. (|T6D due to the Tr '^^'^^HF^y{y)H 
in Eq. ([TT|). The net result for the "hidden" part of x is 



term 



AM V 2d 

This equation is expected to hold for large M. To this should be added the "manifest" part 
given in Eq. (|5p. 

It is important to compare Eq. (^) with the result for x obtained by the exact numerical 
solution for the model based on Eq. (0). This is gotten as an integral over the properly 
normalized radial functions $ (r ),..., \l/2(?") which appear in the P-wave solution of the 
bound state equation [|TU|: 



= A{t)^ r ■ Tpe''\ Ql = -^A{t)^^{r)pe' 



VAtt 



The spinor p labels the grand spin of the bound heavy meson. The choice = +1/2 
corresponds to p = (1, 0)^ The heavy limit bound state wave function in Eq. (|I5|) corresponds 
to the special choice 

<l>(r) oc u{r) , ^i{r) = -$(r) , ^2(?^) = -2$(r) and "^^ir) = . (19) 
The numerical solution to the bound state equations exactly exhibits these relations for 

M,M* ^ oo |ig. 

Equation (|T^ has an interesting li-dependence and vanishes at d = 1/ V2, which actually 
is not too far from the experimental value of this quantity. In Fig. ^ we compare the d- 
dependence of the exact numerical calculation with the perturbative result of Eq. ([171) . It is 
seen that the large M perturbation approach works reasonably well and the gross structure 
of the hyperfine splitting is reproduced. For a detailed comparison of the two treatments it 
is important to note that for fixed M = M* the binding of the heavy meson increases with 
d. In particular this implies that the wave function is only reasonably localized for large 
enough d. As a strong localization is a basic feature of the perturbative approach it is easy 
to understand why this calculation does not yield the exact (numerical) result for small d. 
In fact, as d increases the agreement expectedly improves. However, upon further increase 
of d (at finite M, M*), the numerical solution to the bound state equations shows noticeable 



deviations from the heavy limit relations ([T9|) , which causes the moderate differences at 
larger d. 
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Figure 2: The d dependence of x fo^' ^ = ^* = 30 GeV and d = d! . Solid 
line is the exact numerical calculation. Dashed line is the large M perturbation 
formula given in Eq. ([T^. 

5 Discussion 



We have solved the apparent puzzle associated with the use of a model Lagrangian containing 
ordinary fields for computing the hyperfine splitting parameter x by carefully expanding 
the Lagrangian in powers of 1/M. The key point was the need to preserve the constraint 
■ Q ' = for the heavy vector fluctuation field. 

Of course, such a model Lagrangian (which has been used in many calculations) is not 
exactly QCD. Nevertheless it seems reasonable since it automatically has the correct rela- 
tivistic kinematics and satisfies the heavy spin symmetry at leading order. We have seen 
(Eq. (|ll])) that at next order in 1/M, it predicts the coefficients of many terms which other- 
wise would be unspecified by heavy spin symmetry (even if reparameterization invariance ||T6| 
were taken into account). 

It is amusing to note that these 1/M suppressed terms involve two derivatives and are 
actually more important for the computation of x than the zero derivative term in Eq. (^). 
This is readily understandable since the dynamical scale in this calculation is the binding 
energy, m{B) + m{N) — m{Ab) ~ 620 MeV which is rather large for neglecting light vector 
mesons, higher derivatives etc. [See, for example, Ref. ||17|| .] 



7 



We are regarding the Lagrangian (|^) as an illustrative model rather than as a realistic 
one for comparison with experiment. As indicated earlier it seems necessary to include, 
in addition to finite M corrections, the effects of light vector mesons as well as nucleon 
recoil. The discussion of x this more complicated model and further details of the present 



calculation will be given in a forthcoming publication [15|. 
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